Measurement of the Proton Asymmetry Parameter C in Neutron Beta Decay 
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The proton asymmetry parameter C in neutron decay describes the correlation between neutron 
spin and proton momentum. In this Letter, the first measurement of this quantity is presented. The 
result C — —0.2377(26) agrees with the Standard Model expectation. The coefficient C provides an 
additional parameter for new and improved Standard Model tests. From a differential analysis of 
the same data (assuming the Standard Model), we obtain A = —1.275(16) as ratio of axial-vector 
and vector coupling constant. 
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The free neutron, decaying into electron, proton, and 
anti-electron-ncutrino (in the following called neutrino), 
constitutes a well suited low energy laboratory to study 
the structure of semileptonic weak interactions and to 
probe the underlying symmetries. It is a rather simple 
system without any corrections due to nuclear structure. 
In the V— A description of the Standard Model (SM), it is 
governed by two free parameters only (the Fermi constant 
Gf is considered to be known exactly and CP violation 
can be neglected) : The first element of the quark mixing 
matrix, V u d, and the ratio of axial- vector and vector cou- 
pling constant, A = qa/ gv- On the other hand, there are 
much more experimentally accessible observables, such 
as the lifetime r„ and various correlations between spins 
and momenta of the neutron and the decay products. 
Thus, the problem is overdetermined and permits to test 
the SM description of neutron decay in many ways [l|, Q • 

Prominent examples for angular correlation coefficients 
are the electron asymmetry parameter A and the neu- 
trino asymmetry parameter B, correlating neutron spin 
with electron and neutrino momentum, respectively, and 
the correlation a between the momenta of electron and 
neutrino. Within the SM, each of these parameters al- 
lows to derive A. They have been measured several times 
with increasing precision and arc still in the focus of cur- 
rent research. However, the correlation between neutron 
spin and proton momentum, the proton asymmetry pa- 
rameter C, has not been determined experimentally so 
far. It describes the proton angular distribution Q 

lu(0) d6 = (l + C cos (9) d(9 (1) 

with the angle 9 between neutron spin and proton mo- 
mentum. This Letter reports on the first measurements 
of C. 

A precisely known parameter C permits new cross- 
checks within the SM at low energies as the three decay 
products are kinematically coupled, 



where xc = 0.27484 is a kinematical factor Within 
the SM, C gives an experimentally different access to A 
(gv and gA are assumed to be real here): 



C = x c - 
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C = -x c (A + B) 



(2) 



although with somewhat smaller sensitivity than a or A. 
Physics beyond the SM, e.g. V + A-, scalar- or tensor-like 
interactions, may alter this relation [4, 5]. Thus, C is a 
new input parameter for global analyses of neutron decay 
to constrain the size of possible "new physics" effects. 

A first measurement of the proton asymmetry param- 
eter C using the electron spectrometer PERKEO II has 
been carried out in 2001 [6(. In the following, we report 
on the succeeding experiment of 2004, performed with 
much cleaner systematics and improved statistics. 

PERKEO II was installed at the cold neutron beam 
position PF1B Q of the Institut Lauc-Langcvin (ILL), 
Grenoble. Details on the experimental setup can be 
found in j^, @. The neutron beam was spin polar- 
ized using two super mirror polarizers in X-SM geometry 
[Tot . Polarization P = 0.997(1) and spin-flipper efficiency 
F = 1.000(1) were independent of position and wave- 
length, and stayed constant during the measurement. 

PERKEO II consists of a pair of superconducting coils 
in split pair configuration, generating a magnetic field 
perpendicular to the neutron beam, with a maximum 
of 1.03 T. The neutron spins align parallel to the field, 
that thereby divides the full solid angle into one hemi- 
sphere in and one against spin direction. Each hemi- 
sphere was covered by an electron detector consisting of 
a large area plastic scintillator with photomultiplicr read- 
out. All charged decay products generated in the spec- 
trometer center ("decay volume") were guided onto one 
of the detectors by the magnetic field, realizing full 2 x 27T 
coverage. 

Protons were detected in coincidence with the decay 
electrons alike in refs. 0, |Tl| . The typical energies of elec- 
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tron and proton differ by three orders of magnitude. In 
order to detect both particles with the same detector, the 
protons were accelerated onto a thin (15—30 /ig cm' 2 ) 
carbon foil on negative high voltage (HV). Thereby they 
gained enough energy to release one or several secondary 
electrons from the foil, which were accelerated to the 
plastic scintillator on ground. The decay volume itself 
was kept on ground potential by four grounded grids 
made from AISi wires on each side. Thus all protons 
gained the same acceleration and the assignment of the 
proton to the hemisphere was not affected by the accel- 
erating HV. The detection of the primary decay electron 
was not altered by the HV if the initial kinetic energy 
was above 84 keV, regardless of the emission angle. The 
energy loss of the electrons in the carbon foil was very 
small (~0.5 keV) and constant for all relevant electron 
energies, thus not leading to a systematic effect. The de- 
tection signature was an electron trigger, whose energy 
was recorded, followed by a proton signal. The latter 
arrived delayed due to the finite drift time from the de- 
cay vertex to the region of acceleration. No information 
on the proton energy was available with this method as 
the proton was detected after acceleration and as its drift 
time in the grounded region strongly depended on the po- 
sition of the decay vertex and the emission angle relative 
to the magnetic field lines. 

For each detector, we obtain four spectra, defined by 
particle emission in the hemisphere in (+) or against (— ) 
neutron spin direction: Q ++ (E), Q~ + (E), Q (E), and 
(E). The first sign denotes the electron, the second 
the proton, and E the electron's energy. In order to 
derive the proton asymmetry, we have to integrate out 
the electron (energy and emission direction). With 
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we denote the number of protons that were emitted into 
one hemisphere with the neutron spin pointing towards 
and away from this hemisphere, respectively. In this no- 
tation, the proton asymmetry is 



C 



P - P~ 



P~ 



(5) 



Note that we use only one proton hemisphere and flip 
the spin, thus avoiding systematic effects due to different 
proton efficiencies. The definition ([5|) has an opposite 
sign compared to 0, E3] to maintain the convention that 
a positive asymmetry parameter indicates more particles 
to be emitted in spin direction. 

The finite electron detector threshold, the HV-barrier 
modifying the electron detection at low energies, and 
electron backscattering effects prevent the evaluation of 
the integrals (j4|) in the whole energy range. Further- 
more, the integrals combine events with electron detec- 
tion in two different detectors, without cancellation of 



the related detector efficiencies. Therefore the measured 
Q-spectra were fit by theoretical functions in regions high 
enough to describe the spectra properly and where the 
detector efficiencies for electrons were 1. The spectra 
were then extrapolated to the full energy range and in- 
tegrated, taking into account the finite detector energy 
resolution. 

The theoretical functions to describe the Q-spcctra are 
given in ref [12] with Q 1 ^ = q n F(E) and the Fermi func- 
tion F(E): 
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with 



A and B are the electron and the neutrino asymmetry 
parameter, respectively. P is the neutron beam polariza- 
tion, = v/c. E max = 782 keV is the endpoint of the 
beta spectrum, m c the electron mass. Definition ([6]) is 
separated into two regions (r < 1, r > 1) by the energy 
dependent parameter r, with r = 1 for E = 236 keV. 

The electron detectors were calibrated using six mo- 
noenergetic conversion electron peaks ( 109 Cd, 139 Ce, 
113 Sn, 137 Cs, and two peaks of 207 Bi) covering the rel- 
evant energy range up to 1 McV. They were inserted 
into the spectrometer center regularly. To correct for 
spatial detector effects, the whole decay volume projec- 
tion was scanned two-dimcnsionally several times. Above 
100 keV, detector response was linear within 0.8%. At 
lower energies, the linear relation E = gk + E Q e be- 
tween energy E and analog-to-digital converter (ADC) 
channel k was no longer valid. The detector gain g was 
known with a precision of 0.3% (0.6%) for detector 1 (de- 
tector 2). The energy offsets of the two detectors were 
E oSA = (37.7 ± 1.0) keV and £ off>2 = (40.4 ± 2.1) keV, 
respectively. Their uncertainty is the limiting factor for 
the measurement of C as it directly enters the fit 
functions and thus alters extrapolation and integration 
results. It leads to an uncertainty of 0.82% on C. The 
uncertainties on detector gain and the energy resolution 
of about 15% (FWHM) at 1 McV lead to much lower 
errors (0.38% and 0.12%, respectively). 

A valid event consisted of an electron trigger followed 
by a coincident delayed proton signal. All delayed signals 
within a window of 82 fxs were recorded for both detec- 
tors. The first 40 /is defined the coincidence window. It 
was chosen such that the number of protons with higher 
flight-times was almost negligible (0.03(1)% correction on 
C). Accidental coincidences were measured in the second 
half of the 82 /is window (delayed coincidence method). 
The correction of —0.81(15)% due to accidental coinci- 
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FIG. 1: The spectra to determine C from proton detector 2 were fitted in the regions indicated by the solid lines. The functions 
were extrapolated to cover the whole energy range and finally integrated (filled areas). The one parameter fits describe the 
spectra very well for energies above 240 keV. The agreement between fit and data in the low energy regime can be improved by 
including the detector threshold and a backscattering correction to the theoretical description. We refrain from this approach 
since the uncorrected Q-spectra, as shown in this figure, are needed to determine C. Since the magnetic mirror effect (cf. text 
and acts differently on Q" (electron and proton in same hemisphere) and Q % * , i j (electron and proton in different 
hemispheres), the fit regions are different. 



dences was obtained directly from measured data as the 
timing information for all delayed signals was available 
for all events. 

The fit regions to evaluate the Q spectra were cho- 
sen according to the detailed analysis of the neutrino 
asymmetry parameters [8, 9] generated from same (Q ++ , 
Q~~) and opposite hemisphere (Q H , Q h ) spectra. 
Due to the finite neutron beam width, some particles 
were emitted towards an increasing magnetic field and 
could be reflected (depending on their emission angle). 
For Q ++ and Q~~, this "magnetic mirror" effect in- 
creases for large energies E. Therefore, a smaller fit re- 
gion was used for these spectra. The electron efficiencies 
of the detectors were verified by analyzing the trigger 
functions and found to be 1 above 200 keV. The proton 
asymmetry parameter itself varies slightly (±0.8<7 s t a t) 
with different fit regions and their combinations, lead- 
ing to an additional uncertainty of 0.35%. 

At low energies, there is some HV related background 
left in the spectra (Fig. [T]). In the fit regions, however, 
they can be well described by fits with a normalizing fac- 
tor as the only free parameter, indicating that almost 
no background is left in this region. The remaining un- 
certainty of 0.18% was determined by using the residu- 
als of the fits. Backscattered electrons assigned to the 
wrong detector also distort the spectra below 240 keV. 
At higher energies, backscattering leads to completely 
negligible corrections in the order of 0.01% [l3T |. 

Corrections associated with the design of the spectro- 
meter, "Mirror Effect", "Edge Effect" , and "Grid Effect" , 
are described in They were determined by Monte 
Carlo simulations. For the fit regions of this analysis, 
all associated errors and corrections to C are small. The 



same holds for the influence of the correlation coefficients 
A, B, and a entering via eq. (|6"|). 

Fig. Q] shows the fits and integrals for one detector 
yielding the result 

C = -0.2377(10) stat (24) syst = -0.2377(26). (7) 

A detailed summary of all corrections and uncertainties 
can be found in Table U The second proton detector had 
a much worse proton efficiency due to an inferior carbon 
foil quality, leading to an asymmetry result with am more 
than 3x larger uncertainty. Therefore we only consider 



TABLE I: Correction and errors of the proton asymmetry 
parameter C. The extrapolation uncertainty contributes to 
the statistical and the detector calibration errors. 



Effect [%] 


Corr. 


Err. 


Polarization 


0.30 


0.10 


Flip Efficiency 




0.10 


Data Set: Statistics 




0.44 


Fit Region 




0.35 


Accidental Coincidences 


-0.81 


0.15 


Background 




0.18 


Detector: Gain 




0.38 


Offset 




0.82 


Resolution 




0.12 


Spectrometer: Mirror Effect 


0.01 


0.06 


Edge Effect 


0.26 


0.05 


Grid Effect 


0.08 


0.05 


Correlations: A, B, a 




0.07 


Sum 


-0.16 


1.11 
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FIG. 2: The experimental proton asymmetry C Bxp (E) of pro- 
ton detector 2 and the fit in order to obtain A. Since some 
systematic effects cancel in the asymmetry, the asymmetry 
can be described very well over a large energy range. The fit 
region is indicated by the bold line. 



the precise value here, the second serves as cross check: 
Both values agree within their uncertainties. We also do 
not consider the result of the earlier measurement Q with 
PERKEO II due to its larger systematic uncertainties. 

Our final result for the proton asymmetry parame- 
ter, eq. |(7J), has an uncertainty of only 1.1%. Note 
that it is correlated with the recently published neutrino 
asymmetry parameter B Q as the analyses use partly 
the same data (Q ++ , Q __ ). Our result agrees with 
C*sm = —0.2392(4), the SM expectation calculated using 
eq. © and the world mean value A = -1.2695(29) 0. 
It also fulfills relation ([2]), but the precision of C has to 
be improved for stringent SM cross checks. However, for 
the first time, a proton asymmetry parameter can now 
be included in general parameter space scans of low en- 
ergetic beta-decays to set limits on processes beyond 
the SM. 

We can use eq. §5$ to derive A( int ) = —1.282(21) from 
the integral proton asymmetry C. A more sensitive and 
direct determination of A from the same data can be per- 
formed by a differential analysis of the electron energy E 
dependent proton asymmetry 



Cexp(-E) 



(Q H 



-+) + (Q+- + Q— )' 



(8) 



The dependencies of the Q u on E have been omitted 
here. The fit function is calculated using cq. © and 
replacing the coefficients a, A, and B by their SM ex- 
pressions in terms of A (as given in, e.g., (l2l]). In order 
to avoid systematic effects due to different proton effi- 
ciencies, we again use always the same hemisphere for 
proton detection and flip the spin. Consequently, the 
asymmetry ((5J) contains data from different electron de- 
tectors with different calibration uncertainties. In order 



to estimate the related systematic effects on A, we con- 
sider the properties of the inferior detector. Again, this 
uncertainty dominates the overall error budget. 

On the other hand, systematic effects such as the de- 
tector threshold and the magnetic mirror effect cancel 
in the asymmetry spectrum ((5J) allowing to describe the 
data over a large energy region (Fig. We can extend 
the fit region to up to 700 keV. The lower border is again 
240 keV to avoid systematic effects due to different elec- 
tron detector thresholds or backscattering. The result 



A 



(diff) 



-1.275(6) sta t(15) syst = -1.275(16) (9) 



does not depend on the fit region as results for different 
regions agree within 0.2cr s tat- The overall uncertainty of 
1.2% is dominated by detector calibration (gain 0.68%, 
offset 0.82%) and statistics (0.50%). All other uncertain- 
ties are much smaller. This value agrees with the value 
A.(int) derived from the integral proton asymmetry which 
is an important cross check for the integration and ex- 
trapolation procedure. 

We recommend the value A(diff), eq. ((9]), as result for A 
from our experiment as it is more precise and more direct 
than A(i nt ) . It agrees well with the world mean value A 



14| derived from measurements of the beta asymmetry 
A. It is almost one order of magnitude less precise but 
has been obtained using a new method with different 
systematics: Background is virtually no issue due to the 
coincidence condition. 
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